Quasi-ld quantum helimagnets: The fate of multipolar phases 
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Coupled frustrated spin-1/2 chains in high magnetic fields described within the ferro- antiferro- 
magnetic 3\-J% Heisenberg model are studied by DMRG, hard core boson, and spin wave theory 
approaches. Multipolar phases related to magnon bound states are destroyed (supported) by weak 
antifcrromagnetic (ferromagnetic) interchain couplings J ic . We show that quantum spin nematics 
might be found for LiVCuC>4 whereas for Li(Na)Cu2C>2 it is prevented by a sizeable antiferromag- 
netic J ic . Also for L^ZrCuCU with a small antifcrromagnetic J ic expected triatic or quartic phases 
are unlikely, too. The saturation field is found to be strongly affected even by a relatively small J ic . 



The frustration of magnetic interactions allows for 
new phenomena to emerge. Some spin-chain compounds 
are frustrated magnets in, for instance, the case of fer- 
romagnetic (FM) nearest neighbor (NN) exchange inter- 
action (Ji) competing with an antifcrromagnetic (AFM) 
next-neighbor (NNN) one (J2) [2J. Recent theoretical 
studies indicate that in a high magnetic field, such a J\- 
J2 chain can display e.g. nematic quasi long-range order 
for certain a=— J2/J1 [34a|- This nematic state might 
be thought of as a condensate of two-magnon bound 
states [6j]. Depending on the value of a also three- or 
four- magnon boundstates might condense, resulting in 
a very rich phase diagram with exotic magnetic multi- 
polar (MMP) phases . However, to establish whether 
these MMP phases can be realized in a real material it 
is essential to determine their robustness, i.e. the very 
existence of multimagnon bound states (MBS) with re- 
spect to various interchain (IC) couplings J ic , which can 
be (very) small for certain spin-chain compounds, but 
that is always present, and, as we will show here can be 
detrimental. Such a transition from ID to 2D or 3D can 
be non-trivial. In particular, from quantum mechanics it 
is well-known that bound states are strongly pronounced 
in ID and more rare in 2D or 3D. To understand quanti- 
tatively the role of IC the question arises in which cases 
even a relatively weak IC is still important or even cru- 
cial? Here we address such a problem, namely, the fate 
of MMP states, such as the nematic spin state related to 
MBS at high magnetic fields and consider the spin-1/2 
frustrated isotropic Ji-J2-Heisenberg model mentioned 
above supplemented by various IC [7[. Recently, the ID- 
model and related compounds revealed considerable in- 
terest [2j, |8rl22| . Nowadays it is the mostpopular model 
for edge-shared chain cuprates (see e.g. |2|). Additional 
AFM couplings typically provided by the IC enhance the 
kinetic energy of magnons. Hence, AFM(FM) IC may 
hinder (favor) the formation of low- lying MBS. A FM IC 
may create even new MMP states. An examination of 
real Q-1D systems to predict the changes due to finite 
J ic and to evaluate the chances to detect MMP states in 
certain compounds is of broad interest [a, ll6H2fl l23l - l26| . 



Among edge-shared cuprates LiCu202, the isomor- 
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FIG. 1: (Color) Crystal structure and exchange patterns of 
chain cuprates. The main in- and interchain exchange paths 
are marked by arcs and lines, respectively. Upper: 1/2-unit 
cell projection (right). Planes with green and yellow Cu sites 
are treated as uncoupled. We take into account only per- 
pendicular inplane J ic . Middle: The main 2D IC in the basal 
plane of unshifted chains (left) and with 3D IC (right). Lower: 
Two NN shifted chains in different (ab)-planes. We ignore the 
weak J[ c and keep only the second diagonal J ic [171 ]. 



phic NaCu 2 2 , and LiVCu0 4 =LiCuV0 4 )M (see Fig.Q] 
for comparison the reference systems Li2Cu02 and 
Ca2Y 2 Cu50io with strong IC are shown, too) have been 
proposed to be candidates for quantum-spin nematics, 
i.e. quadrupolar phases derived from 2-MBS. Systems like 
Li 2 ZrCu04 [12|, l22| located closer to the FM-spiral crit- 
ical point might be candidates for the triatic or quartic 
MMP. Besides the strength also the influence of vari- 
ous IC topologies resulting from different arrangements 
of individual chains is of interest. Various types of IC are 
shown in Fig. [TJ The simplest case is given by unshiftcd 
neighboring chains and a predominant perpendicular J^ . 
Here spirals on NN chains are only weakly affected by an 
AFM IC [la] (classically their pitch angle, i.e. the in- 
commensurate (INC) inchain magnetic structure in the 
dipolar phase at ambient external fields remains even un- 
affected). Hereafter, we call this IC the 'unfrustrated IC. 
An effective 2D arrangement of the magnetically active 
Cu 2+ -sites is realized approximately for Li(Na)Cu2Q2, 
where J lc -(0.5 to 1) J 2 ~ 40 to 100 K J2§Hi|. A square 
lattice of unshifted chains, i.e. the 3D case considered in 
Ref. l2J might be realized in LiVCu04, if the IC denoted 
as J3 in Ref. [25| is dominant as compared with the IC in 
the (afr)-plane (see middle of Fig. 1). Another 3D case 
but with shifted adjacent chains is realized for Li2Cu02 
(see Fig. 1) for which H B has been found recently [181 ] . 
J ic «10 K can be extracted from H s (0), from inelas- 
tic neutron scattering studies, and from bandstructure 
results 17 1 . The same order holds also for LioZrCuC^ 
where buckling of the Cu02 chains reduces J ic [14| . 

We used the density-matrix renormalization group 
(DMRG) method [31| with imposing periodic bound- 
ary conditions (PBC) for all directions. In general, 
it is known that this method is much less appropriate 
for D > 1. However, spin systems with up to about 
•s/n x y/n x L = vlO x vlO x 50 sites can be studied 
by taking a proper construction of the lattice block (see 
Ref. UM). We kept m s=s 800 - 4000 density-matrix eigen- 
states in the renormalization procedure. In fact, about 
100 — 300 sweeps are necessary to obtain the GS en- 
ergy within a convergence of 10~ 7 Ji for each m value. 
All calculated quantities were extrapolated to m — > 00 
and the maximum error in the GS energy is estimated as 
AE/Ji — 10~ 4 , while the discarded weight is less than 
1 x 10~ 6 . Under the PBC, a uniform distribution of (Sf) 
may give an indication to examine the accuracy of DMRG 
calculations for spin systems. Typically, (S z ) — S z ot /(nL) 
is less than 1 x 10~ 3 in our calculations. Note that for 
high-spin states [5f t > (nX - 10)/2)] the GS energy 
can be obtained with an accuracy of AE/Ji < 10~ 12 
by carrying out several thousands sweeps even with m w 
100 — 800. Furthermore, we studied the systems with sev- 
eral lengths: L = 16 - 64 (24 - 96) for 3D (2D) and then 
adopting power laws we performed a finite-size-scaling 
analysis to obtain the saturation field h s = gfj, B H s /\Ji\ 
in the thermodynamic limit L — > 00. As a result, we ob- 
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FIG. 2: (Color) Upper (a-c): Saturation field vs IC; lower: 
critical IC for the 3D case with perpendicular IC j lc = 
N ic Jit /\Ji\ vs - inchain frustration a (d,e), where JV ic = 4 is 
the number of nearest interchain neighbors. The pase bound- 
aries are constructed from the 'kinks' in ft a -plots like in (a-c). 
The contour lines in (d) show the values of h s . The 2D case 
with N ic =2 is shown also in (d,e) (dotted lines). 



tain h s with high accuracy. In addition to DMRG we have 
also applied the linear spin wave theory and the hard core 
boson approach [l5j . The latter provides exact results for 
the nematic phase. Some results for the reported magne- 
tization curves have been additionally checked applying 
the exact diagonalization. 

We start with the simplest IC case of parallel chains 
within a plane (2D) or within a square-lattice arrange- 
ment of chains (3D) and a single perpendicular IC 
J^ =j'rr I J\\ (see Figs. 1,2). A typical magnetization curve 
for the nematic phase is shown in Fig. 3. The height of 
the magnetization steps AS Z = 2, only, is its signature 
whereas in the 1-magnon phase AS Z = 1 holds. Note 
that a rather weak critical IC of a few percent will re- 
move the nematic phase in favor of the usual field in- 
duced "conic" phase. The saturation field h s of the INC 
phase on the 1-magnon side is described exactly already 
within spin wave theory: 

h s = 9^Hj\,h\ = 2a- 1 + 0.125/Q- + j ic , (1) 

where j ic =N ic j ic =2{D - l)j ic for D=2,3. Using the Curie- 
Weiss temperature 9 C „ = -(Ji + J 2 + 0.5A ic J ic )/2 
which determines the high-temperature spin susceptibil- 
ity x(T) - 1/(T - e cw ), the IC in Eq. (1) can be elim- 
inated. Then in this 1-magnon phase one has a simple 
expression to extract J\ or a from H s and O cw : 

3MbH s (J 1c ) + 49 CW (J 1C ) =G = |Ji| [1 + l/(8a)] . (2) 



A relation G = \Ji\f(a) is valid also in other regimes, 
where / is a simple function affected by the type of the 
IC and the region of the MMP phase diagram. In the 
1-magnon regimes valid for LizCuC^ and Ca2Y2Cii50io 
(see below) we have / = 2(1 — a) [18(. For the nematic 
phase in full accord with the DMRG we obtained numer- 
ically exact values of h B using the hard-core boson ap- 
proach [15|- Expanding the 2-particle Green's function 
in powers of j l0 <C 1 we arrive at analytical expressions: 

h s (a) a hf(a) + j\./2 + N icV (a)j? c + O (j, 4 c ) , (3) 

hl°(a) = 2a -1 + 0.5/(1 + a), (4) 

where 77(a) = (1 + a)(3a 2 + 3a + 1)/[2(1 + 2a) 2 ] « 
5/3 + 3a/2 . The expansion coefficients for higher or- 
der terms differ for D=2 and 3 reflecting subtle natural 
differences in the quantum fluctuations. A detailed dis- 
cussion will be given elsewhere. Comparing the Eqs. (1) 
and (3) we stress the presence of nonlinear IC terms and 
two times smaller linear term in the nematic phase as 
compared with the usual one-magnon phase. The in- 
spection of Figs. 2(d,e) shows that the critical IC of the 
nematic region reaches a maximum as a function of a. 
Its position a « 1.103 and height ji c /4 ~ 0.071688 are 



close to the feature shown in Fig. 6 of Ref. [24 . 

Keeping the linear IC term in Eq. ([3]) we arrive at 



[12j. The small ID- value of H, w 4.2 to 5.9 T, only, 
clearly shows the importance of J ic in the vicinity of the 
critical point where H s — > 0. The validity of Eq. @ is 
guaranteed by j ic = 0.0264756 (i.e. 7.2 K) to 0.023659 
(i.e. 6.5 K) exceeding well j ic ,„=0.009539 to 0.013458 de- 
duced from the 2D-phase diagram. The j iccr would allow 
H s = 6.2 T to 8.9 T in the triatic phase at most. The esti- 
mated weak IC given above is close to L(S)DA+£/ results 
14j . Thus, we are clearly outside the triatic or quartic 
phases and deep enough in the usual dipolar phase. The 
1-magnon picture holds also for Li(Na)Cu202 where j iccr 
for the nematic phase is exceeded by a factor of four. 
In fact, from the LDA derived J ic and J\ we estimate 
j ic « 0.7 well above 0.170 j lr ..„ = 0.1704 taken from Fig. 
2 for a sa 1 for LiCu9Q9[28l|. Similarly, using the empir- 
ical values for NaCu202: Q cw = —41 K taken from the 
l/x(T) data at 350 < T < 400 K, a = 1.9, J x « -48 K, 
and g a = 2.06, we predict H s ~ 155 T. From Eq. (1) we 
estimate j ic ~ 0.73 well above jL„ ra 0.172 (see Fig. 2) 
and again close to LDA-results [30j. 

Now we turn to the case of diagonal IC as in Li2Cu02 
(see Fig. 1). For a strong enough IC the MBS are su- 
pressed and only 1-magnon excitations survive at high 
fields (see Fig. 4(b)). Then h s reads 

K = gi* B H./\j 1 \=j lB + j , Ul if j ic > %„ h = N iJic (6) 



gHvH s + 20c 



K = | Jx | [a + 0.5/(1 + a)] 



(5) 
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Thus, G and K composed from IC affected quantities de- 
pends itself on single chain properties, only. In most cases 
it is easier to determine the latter theoretically whereas 
G or K can be found from experiment. 

Applying Eq. ^ to Li2ZrCu04 we predict C w = 
93.3 K for a preliminary experimental value of H s of 
about 13 T J3i| using J x = 273 K, a = 0.29 to 0.3 
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FIG. 3: (Color) Magnetization vs. applied external field for a 
2D arrangement of chains with a direct AFM IC as modelled 
by four chains with N = 24 sites in each chain for different 
IC ji C = Ji /|Ji| and an inchain frustration a — 1/2. Note 
that each point shown represents a step. The critical value of 
j ic /2 amounts 0.088 for D=3 and 0.094 for £>=2. 
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FIG. 4: (Color) 3D phase diagrams for diagonal IC. h s (given 
by contour lines) as a function of the normalized IC j ic and 
the inchain frustration a. Lower(Upper): (un)shifted chains 
with ±36/2 (±b) IC, where b is the inchain Cu-Cu distance. 



where iV ic = 8. Eq. © is valid for j lc {a) > (4a - l)/9 = 
j" t for j' ic = 0. We obtained from INS data [l7J|: 
Ji C w 9K > J^ (0.332) = -0.0364Ji «8.2 K. Then 
H H depends solely on the IC which can be directly read 
off from H S (T = 0) 01 . For weaker IC /i s depends 
also somewhat on a. Our results also suggest that in 
this intermediate INC-phase above a second critical IC 
Ji*2 ~ 0.0109 only specific INC 1-magnon low-energy ex- 
citations exist, see Fig. 4. Below j" 2 3-MBS are recovered 
as low-energy excitations and only this narrow region is 
ID-like. The commensurate (C) phase behaves like an 
ordinary 3D antifcrromagnct despite its seemingly quasi- 
1D nature. For unshifted NN chains (see Fig. 1), the 
C-phase is missing. There is only one critical IC separat- 
ing INC 1-magnon excitations from 3-MBS. 

Generally, 9 = — 2@ cw / J\ provides a useful constraint 



l = a + l3 + 9 + (D- l)\j i0t± +2j^ j h 

f 



(7) 



for the exchange integrals, where in addition to a further 
inchain and IC quantities j3 = — J3/J1, j le>± = —J lc ,x/Ji) 
and j iCil = —Jic.t/Ji for various diagonal couplings have 
been introduced. From Eq. ([7]) a stringent constraint for 
a for LiVCu0 4 follows. From l/x(T) at 500 K < T < 
650 K a small FM value ©cw = 7.4 K can be fitted from 
the data shown in Refs. |25j, |33J. The FM ordering of spi- 
rals observed by neutron diffraction can maintained both 
by a FM inplane and a perpendicular AFM 3D IC (see 
Fig. 1) Then we may assume for the sake of simplicity 
that both IC almost cancel in Eq. ([7]) and ignore the small 
terms j3, and j i0 in zeroth approximation. The weakness 
of the total IC is suggested also by the weak magnetic 
moment m = 0.3/x B at T= [34j pointing to strong quan- 
tum fluctuations. A first estimate [33J of the dispersion 
of the INS peaks as well as of the INS intensity above 
10 mcV yields a reasonable description for a=0.8 and 
J\ =— 73.2 K and J3—— 0.07 K as suggested by a mapping 
of a five-band Hubbard model (exact diagonalizations) 
for a CugOi2 cluster onto a six-site J1-J2- J3-Heisenberg 
ring. Hence, we conclude, that for a w 0.85 ± 0.1 and 
J 1 » -80 K nematics could be observed in LiVCuGi. 
Using Eq. (J3j) we see that this is also in reasonable ac- 
cord with the experimental -ff s -data 20j,[25| since our ID 
set gives for the above mentioned rough estimate about 
48 T for g — 2 and H directed along the hard axes and 
about 42 T for the weak axis with g = 2.3. Note that 
our J's differ strongly from those given inappropriately 
in Refs. |25|, [26j. The latter have been employed unfor- 
tunately in Refs. IfJ l20j as 'empirical' input parameters. 
More detailed studies including explicitly the IC are nec- 
essary to refine all J's. In particular, the quantum effect 
of J ic on the pitch angle is of interest. A more detailed 
study of FM IC as well as of the J1-J2-J3 inchain model 
will be considered elsewhere. A challenging point is also 
to find MP phases beyond the 2-magnon Bosc condensa- 
tion [6j triggered by a FM IC. 



To summarize, the crucial role of realistic AFM IC in 
Q-1D hclimagnets has been demonstrated. The rich and 
exotic physics of multipolar phases recently predicted for 
single chains is very sensitive to the strength of the IC. It 
can be readily eliminated by a weak AFM IC especially 
for triatic, quartic, etc. phases. For most Cu02 chain 
systems studied so far, except probably LiVCu04, where 
a nematic phase or some other exotic phase might be ex- 
pected, the AFM IC is too strong to allow for multipolar 
phases. The examination of various weak FM IC as well 
as of anisotropy effects in real materials is under study. 
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Note that the ordered moment for LiVCuCU is three 
times smaller than that for Li2CuC>2 or Ca2Y2CusOio. 



